(N : Phenomenology of SIDIS unpolarized 

3 '• cross sections and azimuthal asymmetriesli] 






09 ' Vincenzo Bar one 

(N 



^ 



(N 



Di.S.I.T., Universita del Piemonte Orientale "A. Avogadro' 



P^, INFN, Gruppo Collegato di Alessandria, 15121 Alessandria, Italy 

hi- 



> 

!~! ! Abstract 

\q 

^^ ' I review the phenomenology of unpolarized cross sections and azimuthal asymme- 

^D ■ tries in semi-inclusive deeply inelastic scattering (SIDIS). The general theoretical 

framework is presented and the validity of the Gaussian model is discussed. A 

brief account of the existing analyses is provided. 
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c3 : 1 Introduction 



It has been known since the early years of quantum chromodynamics that az- 
imuthal asymmetries in unpolarized processes are generated by gluon radia- 
tion and splitting. The observation of these asymmetries was in fact proposed 
as a test of perturbative QCD \T\. The first experimental results date back 
to the Eighties and showed that unpolarized semi-inclusive deeply inelastic 
scattering (SIDIS) indeed displays non vanishing azimuthal modulations of 
the type cos (j)h and cos 20/^ [21 [3] . 
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Figure 1: Transverse kinematics in the photon-nucleon coUinear frame. 

While the pQCD processes account for these asimmetries at large trans- 
verse momenta, in the low-Ph± region it is the intrinsic transverse motion of 
quarks that plays a key role (for a recent review see [4j). There are two non- 
perturbative sources of azimuthal asymmetries in unpolarized processes: 1) 
the Cahn effect [HI E] , a purely kinematic correction due to the quark trans- 
verse momentum: 2) the Boer- Mulders effect p], arising from the correlation 
between the transverse momentum and the transverse spin of quarks inside 
an unpolarized nucleon. 

2 Theoretical framework 

We will consider unpolarized SIDIS, /(£) + N{P) -^ Z(f ) + h{Ph) + X{Px), 
in the current fragmentation region. At low Ph± the cross section of this 
process factorizes into a transverse-momentum dependent quark distribution 
(TMD) and a fragmentation function [8]. 

Let us consider the transverse kinematics in a reference frame where the 
virtual photon and the nucleon are coUinear (we call it 7*A^ collinear frame, 
see fig. [l]). The transverse momenta are denoted as follows: k±, transverse 
momentum of the initial quark; p±, transverse momentum of the hadron 
w.r.t. to the fragmenting quark; Ph±, transverse momentum of the hadron 
w.r.t. to the 7*A^ axis. Neglecting l/Q"^ corrections, these momenta are 
related by 

P±^ Ph±- Zhk±. (1) 

The SIDIS structure functions are actually expressed in terms of the quark 
transverse momentum fc-r defined in the frame where the produced hadron h 
and the nucleon are collinear {hN collinear frame). However, if one neglects 
1/Q^ corrections, the transverse momenta in the 7*iV and hN collinear frame 



coincide: k± c:^ kx- 

The general expression of the cross section for unpolarized SIDIS is [U [TU] 
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{(1 - y + 2?/') Puu,T + (1 - y) Fuu,L 
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+ (2 - y)^l^y cos0,F™J*'^ + (1 - Z/) cos20,F™j2^''} . 



The azimuthal asymmetries are defined as: 
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and similarly for A'^"^'^'^'^. 

In the extended parton model, the structure functions appearing in eq. ([2]) 
are given by convolutions of transverse momentum dependent quark distri- 
butions and fragmentation functions. Up to order 1/Q one has 
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where C[w f D] is a transverse momentum convolution. 

Equation (|5]) represents the Cahn contribution to the cos0ft structure 
function. It originates from the elementary lepton-quark cross section in 
presence of transverse momenta, which reads 
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Equations ([6]) and d?]) are the Boer-Mulders (BM) contributions to the 
azimuthal modulations of unpolarized SIDIS. They couple the Boer-Mulders 
function hi to the Collins function Hi [llj, which describes the fragmen- 
tation of transversely polarized quarks into an unpolarized hadron. Notice 



that the BM effect contributes to cos0/i at order 1/Q {i.e., at twist three) 
and to cos 20/1 at leading twist. 

At twist three one should also take into account quark-gluon interactions, 
which give rise to the so-called "tilde" distributions. Thus F^^'' gets the 
additional term [10] 
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We recall however that TMD factorization is not proven beyond leading twist, 
so eq. OH]) must be taken with a grain of salt. 

The Cahn effect produces at order l/Q"^ (twist four) a further contribution 
to the cos 20/j structure function. 
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This expression incorporates only part of the l/Q"^ kinematic corrections. 
Besides them, there are also dynamical l/Q^ corrections arising from twist- 
four TMD's and fragmentation functions. Therefore, -^tJx/cahn i^ J^st an 
approximate estimate of the full (9(1/Q^) structure function. 

In summary, ignoring unknown higher-twist terms, the two azimuthal 
asymmetries of unpolarized SIDIS are symbolically given by: 

(COS0,,) = -Cahn + -BM, (cos 20/,) = BM + — - Cahn. (11) 

Q Q Q^ 

3 The Gaussian approach 

In most phenomenological analyses, the TMD's and fragmentation functions 
are parametrized by Gaussians: 

/i(a;,fci) =^/i(a;)e-^i/^ D^{z,p\) = BD^{z)^-^'^^^ . (12) 

Here x is the fraction of the nucleon light-cone momentum carried by the 
quark, and z is the fraction of the light-cone momentum of the struck quark 



carried by the final hadron. These two variables coincide with xb and Zh, 
respectively, modulo l/Q"^ corrections. We stress that there is a strong as- 
sumption - with no fundamental basis - behind the Gaussian model, namely 
that distributions factorize in x and k±. 

In principle, one has to distinguish the Gaussian widths k']_ and p^ from 
the average squared momenta, defined as 

{kl)^ld'k^klMx,kl), {pl)^JdWplD^{z,pl). (13) 

Average values coincide with the Gaussian widths, {k\) = fcj_, {p\) = pj_, 
only if we integrate over k± and p± between and cxd. 

Bounds on X = k~/P~ and k± can be obtained in the quark parton 
model by inserting a set of physical intermediate states in the quark-nucleon 
amplitudes [l2] . Since the momentum of these states is 

the condition M^ > yields 

kl<x{l-x)M^. (15) 

This is the bound on the truly intrinsic quark transverse momentum [13] (for 
recent discussions see [HJIIS]). Numerically the upper limit ( TT5|) is k^ < 0.25 
GeV^. Clearly, the average squared momentum (k^) must be smaller than 
this value. On the other hand, various phenomenological analyses give much 
larger figures, (kl) ~ 0.25 - 0.40 GeV^. 

Is there a contradiction between these results? The answer is that the 
bound ( IT5l) actually refers to a "static" nucleon {Q^ = 0), whereas the value of 
(k^) extracted from experiments effectively accounts for the "non-intrinsic" 
transverse momentum, radiatively generated at a given Q^. As a conse- 
quence, {k\) is not a fixed, universal, quantity. 

Consider now the unpolarized and angle-independent structure function: 

Fuu = E e' / d'fc± / d'p± S\p^ + z^k^ - Ph±) f^{xB, kl) Dl{zh, pi) 

(16) 
Using Gaussian parametrizations and integrating between and oo, one gets 



Fuu = T.^l fii^B) Dlizn) """^^^ (17) 
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Figure 2: The ratio R{Ph±) = do'uu{Ph±)/(io'uu{0): Gaussian model vs. 
CLAS data (from p]). 



with 



Pl^ = pl + 4kl (18) 

This is a popular relation appearing in many analyses (where it is often 
quoted as a relation between average momenta). It must used with some 
caution. In fact: i) it is invalidated by a possible cutoff on k"^; ii) it has 
(complicated) l/Q^ corrections, in) due to experimental cuts, P\j_ differs 
from the measured {P\x)- 

4 Transverse- momentum dependence of cross 
sections 

The Gaussian Ansatz for the unpolarized SIDIS cross section, 
d(T(7t/(Ph±) da[/c/(0) .pij-^ 



(19) 

h± ^^ ^-^ h± 



dzdPl, dzdPl 



has been tested by Schweitzer, Teckentrup and Metz [16] in the CLAS kine- 
matics {xb = 0.24, Zh = 0.34, g2 = 2.4 GeV^) [17]. As shown in fig. El the 
description of data is very good, with a Gaussian width P|_|_ = 0.17 GeV^. 

The results of a Gaussian analysis of the Plx distributions measured by 
COMPASS are displayed in fig. [3] (left) [18]. The fitted P|7 is found to have 
a mild dependence on xb and to increase with Q^. A check of the relation 
f lT8|) has been also performed, with the result shown in fig. [3] (right). One 
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Figure 3: Left: P^^ distributions fitted by Gaussians (COMPASS data) [18 
Right: Fitted ^ (called (p^) in the figure) vs. zl (COMPASS data) [H].' 
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Figure 4: Left: Dependence of {Ph±) on the lepton-nucleon cm. energy 
squared of various experiments (from [I6]). Right: Dependence of (-P|x) on 
the photon-nucleon cm. energy squared W^ (COMPASS data, from [T8]). 



sees a clear departure from a linear dependence of P^ 
and dot-dashed curves in the figure are fits of the type 



f^^ on zf^. The dashed 
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Finally, it has been noticed in [16] that the average transverse momentum 
of hadrons extracted in different experiments (JLab, HERMES, COMPASS) 
has an approximately linear rise in s (fig. HI left). This energy dependence is 
confirmed by a COMPASS analysis [18], which shows that {Pl±) increases 
with W"^ (fig. m right), which may be an indication of a possible broadening 
of (fci) and (pi). 



COMPASS Deuleran 




Cahn 

Boer-Mulders 

Cahn+BM 

0.3 0.4 0.5 0.6 0.7 0.1 0.3 0.5 0.7 0.9 

z PT(GeV) 



Figure 5: The fit of |2T| to tlie COMPASS preliminary data on {cos2(f)h) witli 
a deuteron target. 



5 Azimuthal asymmetries 



Tlie first analysis of the cos (ph asymmetry in SIDIS in terms of the Cahn effect 
was presented in 2005 by Anselmino et al. [19]. They fitted the EMC data 
[2] with a Gaussian model, obtaining the values fc^ = 0.25 GeV^, p\ = 0.20 
GeV^, but did not consider the Boer-Mulders term, which is of the same 
order as the Cahn contribution. 

The cos 20/1 asymmetry was extensively studied in [201 EI]- Both the 
leading-twist Boer-Mulders contribution and the twist-4 Cahn term were 
included in the analysis. Signs and magnitudes of the Boer-Mulders function, 
extracted from a fit to HERMES |22] and COMPASS [23] preliminary data, 
were found to be in agreement with the theoretical expectations (based on 
the impact-parameter approach, lattice studies and large Nc arguments): 
hi- ~ 2A-L.", hi'' ~ -/iV. 

In particular, it was pointed out that, since the Cahn effect is the same 
for TT"'" and 7r~, and the favorite and unfavorite Collins functions are related 
by Hi^''" ^ -i7^^"°^ as shown by the fits to the Collins effect in SIDIS, a 
signature of the Boer-Mulders effect is (cos20/i).n- > (cos20/i)^+. Another 
interesting output of the analysis in [201 121] is that the Cahn contribution 
turns out to be relatively large in spite of being (^(l/Q^). An example of the 
fits presented in [21] is shown in fig [51 

The HERMES and COMPASS Collaborations have recently presented 
new preliminary data on azimuthal asymmetries in unpolarized SIDIS |24[ 



While the two measurements are in fairly good agreement for (cos0/j), 
the COMPASS result for the cos 20/j asymmetry is systematically larger than 
the corresponding HERMES result. 

The CLAS Collaboration at JLab has released final results on the az- 
imuthal modulations [17]. A large discrepancy is found between the data 
on {cos(j)h) and the Cahn effect predictions (extrapolated from other anal- 
yses). The (cos20/i) signal is non zero only in the small- z/j region, where 
target fragmentation effects are important and may affect the interpretation 
of data. 

From the theoretical viewpoint, the situation is quite unclear, as only 
partial analyses are presently available. If we simply extrapolate them, we 
find that the Cahn contribution to (cos0/j) is huge and largely overshoots 
the data, and that the Boer-Mulders contribution to (cos2(/)/i) is small and 
apparently goes in the wrong direction, in the sense that it has an opposite 
sign in {cos (f)h) and {cos2(f)h), and therefore should produce a {coscph)^- 
smaller than {cos (f)h)iT+, whereas the data show a similar tt^-tt^ pattern for 
both asymmetries. 

We have already said that the intrinsic transverse momentum of quarks is 
kinematically bounded. We expect that the radiatively generated transverse 
momentum should be limited as well. The effect of an upper bound on 
k± has been recently explored by Boglione, Melis and Prokudin [26]. They 
derived a cutoff of the form fc^ < r]{xB)Q'^ in the 7*A^ center-of-mass frame. 
The consequence of setting this bound is that the Cahn contribution gets 
largely suppressed, while the BM contribution remains almost unaffected. 
This correction goes in the right direction, but further work is needed to 
establish a frame- independent condition on k±. 

6 Conclusions 

There is nothing magic about the Gaussian approach. It is just a parametriza- 
tion, which happens to work rather well for cross sections at low Ph±- Since 
this model lacks a solid basis, one should not be surprised if some simple 
regularities based on it turn out to fail. 

The parameters k']_ and p^ in the Gaussian distributions are likely to be 
Q^ (and W'^) dependent. This means that they cannot be fixed once forever 
and used anywhere. One should do as for normal PDF's: take a TMD at 
small Q^ from some model, evolve it and fit the data. 



Although we have clear signals of relatively large azimuthal asymmetries, 
the situation is still unsettled. A combined state-of-the-art fit of (cos0/i) and 
(cos 20ft) would be highly desirable in order to understand the origin of these 
observables. We should be prepared to the possibihty that the scheme (ITTl) 
might not work. We recall that the corrections to this scheme are of three 
types: 1) "genuine" [i.e. "tilde") twist-3 contributions (originating from 
quark-gluon correlations); 2) further 1/Q^ kinematic terms; 3) dynamical 
twist-4 contributions. 

A lesson we have learned is that higher twists are relevant in the presently 
explored region of SIDIS. Clearly, it would be important to disentangle them 
from leading twist contributions. In order to do so, a wider Q^ range is 
needed - an important task for future experiments. In the meanwhile, the 
l/Q^ contributions should be fully worked out, computing both the kinematic 
corrections and the target mass effects. 
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